Bragg spectroscopy of a strongly interacting Fermi gas 
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We present a comprehensive study of the Bose- Einstein condensate to Bardeen-Cooper-Schrieffer 
(BEC-BCS) crossover in fermionic ^Li using Bragg spectroscopy. A smooth transition from molec- 
ular to atomic spectra is observed with a clear signature of pairing at and above unitarity. These 
spectra probe the dynamic and static structure factors of the gas and provide a direct link to two- 
body correlations. We have characterised these correlations and measured their density dependence 
across the broad Feshbach resonance at 834 G. 

PACS numbers: 03.75.Ss,05.30.Fk,03.75.Hh 



Strongly interacting ultracold Fermi gases provide new 
opportunities for the study of pairing and superfluid- 
ity. Magnetic field Feshbach resonances allow precise 
tuning of the interactions between fermions in differ- 
ent spin states and this has led to the realisation of 
long lived molecular Bose-Einstein Condensates (BECs) 
P, H, S 0, S| and fermionic condensates [1, [3|- Since 
their first realisation, the properties of these gases in the 
BEC-BCS (Bardeen-Cooper-Schrieffer) crossover have 
attracted much attention [1, i, [l^, lili III III, Q III . Re- 
cent experiments using radio frequency (rf) spectroscopy 
have been been particularly fruitful yielding the pair size 
16 1 , low energy excitation spectrum and pairing gap 
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One tool that has not yet been applied to these systems 
is Bragg spectroscopy. This differs from rf spectroscopy 
in two important ways. Firstly, Bragg scattering does 
not change the internal atomic states so there are no final 
state interaction shifts. Secondly, Bragg scattering can 
transfer momentum larger than the Fermi momentum, 
hkp, whereas the momentum transferred by an rf photon 
is neghgible (<C tikp). As such Bragg spectroscopy can 
probe a wide region of the excitation spectrum. 

Several theoretical studies of Bragg scattering in Fermi 
gases have predicted si gna tures related to pairing and su- 
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perfluidity 
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22( 1 . Bragg spectroscopy has pre- 
viously been used to characterise the dynamic and static 
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structure factors of bosonic condensates [231 . 
quantities are linked to the two-body correlation func- 
tion g^'^\r) via the Fourier transform [2^ [2^. In an 
interacting two component Fermi gas, strong pair corre- 
lations may exist and these can be probed with Bragg 
spectroscopy. Recent Bragg experiments have studied 
boson/boson molecules near a Feshbach resonance 12711 . 
Ramsey interferometry with noninteracting fermions [28| 
and molecular condensation [2§|, but no work has yet 
investigated correlations in the BEC-BCS crossover. 

In this letter, we report Bragg scattering experiments 
around the broad s-wave Feshbach resonance at 834 G in 
an ultracold gas of ^Li. Bragg spectra reflect the compo- 



below the Feshbach resonance, pairs and free fermionic 
atoms near unitarity, and free fermions above resonance. 
Our spectra show features present in the dynamic struc- 
ture factors calculated by Combescot et al. [2l[. Strong 
pair correlations are observed at unitarity which decay as 
the density is lowered highlighting the many body nature 
of pairing on the BCS side of a Feshbach resonance. 

Our experiments use a gas of ^Li prepared in an equal 
mixture of the \F = 1/2, = +1/2) and |l/2,-l/2> 
states evaporatively cooled in a single beam optical dipole 
trap [30] . This yields a highly degenerate cloud with N — 
1.3 X 10^ atoms in each state at a magnetic field of 835 
G. The field is then ramped in 100 ms to a desired value 
and held there for both Bragg scattering and imaging. 
The imaging light is tuned to be resonant with atoms in 
the |l/2, 1/2) state at each magnetic field. On the BEC 
side of the Feshbach resonance, molecules can be imaged 
using the same frequency light as atoms as their binding 
energy is much lower than the transition linewidth. 

All experiments presented here were performed at the 
lowest temperature achievable on our experiment. At 
unitarity we quantify this by the empirical temperature 
T obtained from one-dimensional fits of an ideal Fermi 
gas profile to our measured density distributions inte- 
grated over two dimensions T provides an estimate 
of the true temperature via T ~ T/ {Tpy/l + (3), where 
Tp = Ep/kB is the Fermi temperature and /3 = —0.58 
is a universal constant 0, [3l[. We find T = 0.1 and 
T/Tp ~ 0.07. Temperatures of clouds at other magnetic 
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fields may be calculated from the entropy 
on the BEC side of the resonance the temperature can be 
estimated from bimodal fits to expanded molecular den- 
sity profiles. At 700 G our measured condensate fraction 
of ~ 90% is consistent with T/Tp < 0.1 at unitarity, 
following an isentropic magnetic field sweep 32 1. 

Bragg scattering is achieved using two counter- 
propagating laser beams with a small tunable frequency 
difference S to create a moving periodic potential. A 
two photon Bragg event transfers momentum hq = 2hkL 
to the scattered particle. Ignoring Doppler and inter- 
sition of the gas, being dominated by bosonic molecules action shifts, the Bragg resonance frequency is given by 
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FIG. 1: (a) Schematic of Bragg scattering experiments. Bragg 
scattering of molecules from (b) a pre-expanded molecular 
BEG and (c) a trapped molecular BEG, both at 750 G. The 
field of view for images (b) and (c) is 650 fim by 485 /im. 
Elastic collisions between scattered and unscattered particles 
distort the resultant cloud in the trapped case. 



S = 2hk\/m, where /cl is the wavevector of the laser 
and m is the mass of the scattered particle. In a Fermi 
gas near a Feshbach resonance, these particles can be 
atoms (with mass M), tightly bound molecules or cor- 
related pairs (mass 2M). Bragg scattering therefore 
allows one to distinguish molecules/pairs from atoms. 
The relative fractions of each depends on the temper- 
ature, T/Tp, and interaction parameter, 1/kpa, where 

Lo is the geometric mean trap- 
ping frequency and a is the s-wave scattering length. For 
^Li the Bragg resonant frequencies are 52m/'2,tt — 147 
kHz for molecules and 5m /^tt = 294 kHz for atoms. 

The Bragg lasers are applied near perpendicular to the 
long axis of the dipole trap. Fig. [T] (a), and are lin- 
early polarised parallel to the magnetic field direction. 
To demonstrate free particle scattering, Fig. [1] (b) shows 
an absorption image of molecules scattered from a molec- 
ular BEC {5/2tt = 140 kHz) at 750 G which has expanded 
for 3 ms before the Bragg pulse. The pulse duration tbt^ 
and intensity were chosen to yield a tt pulse for resonant 
molecules. A thin slice of the cloud has been removed 
demonstrating the high momentum selectivity of Bragg 
scattering. With longer Bragg pulses we have observed 
several Rabi cycles of scattered molecules (and atoms). 

Bragg scattering of pre-expanded low density clouds 
however, cannot reveal the effects of interactions, so we 
now focus on scattering of trapped gases. This can be 
complicated by elastic collisions between scattered and 
unscattered particles 



34l |. These result in distorted and 



asymmetric atomic distributions in which it is generally 
not possible to discern a spatially separated scattered 
cloud. For Bragg spectroscopy of trapped gases, we apply 
a short Bragg pulse {tbt — 40 /is) and switch off the 
trapping potential immediately following the pulse. Both 
of these events happen on a timescale short compared 
to the inverse trapping frequencies {fr = 320 Hz, = 
24 Hz). The cloud then expands for 4 ms before being 
imaged. Fig. [T] (c) shows the outcome of this sequence 
for a trapped molecular BEC at 750 G. 



To analyse such images, we quantify the effect of the 
Bragg pulse using the centre of mass displacement X. 
The momentum imparted to the cloud P(g, 5) by a Bragg 
pulse with wavevector q and energy M determines the 
resultant centre of mass velocity '^^^^''^^ = P{q,S)/m. 
For q > kp and times short compared to th e inverse 
trapping frequencies this is given by [26 
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dt 
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where flgr = (r'^/4A)^//i/2//s is the two-photon Rabi 
frequency, Ii and I2 are the intensities of the two Bragg 
beams. Is is the saturation intensity, F is the natural 
linewidth of the transition and A is the Bragg laser de- 
tuning. A/27r ~ 1.5 GHz is chosen to be large com- 
pared to the 80 MHz splitting between the |l/2, 1/2) and 
1 1/2, —1/2) states so the coupling to each state differs by 
only 5%. S{q,6) is the dynamic structure factor which 
describes the response of the gas to excitations, fi^^ 
is twice as large for molecules/pairs as it is for atoms 
(molecules have twice the polarisability) however, this 
is exactly compensated for in ([T]) by the factor of two 
mass difference. Thus X{q,6) has the same factor be- 
fore the integral for both atomic and molecular (paired) 
gases. Additionally, elastic collisions preserve centre of 
mass motion, so X{q,5) is a good measure of the Bragg 
signal throughout the BEC-BCS crossover, independent 
of cloud shape and scattering length. 

We evaluate X{q, 6) by integrating images such as Fig. 
[T] (c) horizontally and finding the centre of mass of the re- 
sultant line profile. This is taken relative to the centre of 
mass of a reference image obtained with no Bragg pulse. 
Fig. [5] shows Bragg spectra obtained in this way for a 
range of magnetic fields on either side of the Feshbach 
resonance. At 750 G and 780 G {l/kpa = 1.4 and 0.8, 
respectively) the spectra are dominated by a large peak 
near 150 kHz corresponding to scattering of molecules 
from a molecular condensate. In this regime the molecule 
size is smaller than the inverse excitation momentum, 
a < 1/q, so molecules behave like single particles of mass 
2AI. The height of this peak is indicative of a large con- 
densate fraction. The spectra are slightly asymmetric 
showing some excitation of free atoms. Closer to reso- 
nance at 820 G {1/kpa = 0.2) the molecule/pair peak 
decreases as a > 1/q and the spectra merge with the 
continuum of free atomic excitations. At unitarity (835 
G) the spectrum has significant components due to both 
pair and atom scattering. This provides clear evidence of 
a substantial paired fraction in a trapped unitary gas. At 
higher magnetic fields (850 G, 1/kpa = —0.2) the pair 
signal remains strong. For 890 G and 991 G, 1/kpa = - 
0.5 and -1.2 respectively, the pairing feature drops off and 
the spectra approach that of an ideal Fermi gas. These 
spectra show the smooth transition in S{q, S) from being 
dominated by molecular excitations below the Feshbach 
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FIG. 2: Bragg spectra showing X{q,S) for trapped Fermi 
gases across the BEC-BCS crossover. Magnetic field and 
(l/kpa) for each spectrum are given in the legend. The inset 
shows the 750 G and 991 G spectra along with the calculated 
X{q,S) for an ideal Fermi gas and molecular BEC at 750 G. 



resonance to atomic excitations above resonance, in good 
agreement with the predictions of ref . [2l[ . 

The inset of Fig. [U shows the full 750 G and 991 G 
spectra along with the calculated X{q,5) for a T = 
ideal Fermi gas and Thomas Fermi molecular BEC {amoi 
= 110 nm at 750 G). These limiting cases are found using 
the impulse approximation for S{q,d), valid for large q 
37| . The bosonic molecular condensate response is much 
narrower and higher peaked than the Fermi gas response. 
The relative heights of the calculated responses are set 
by requiring the area under the bosonic curve to be twice 
that of the fermionic curve 21]. Both the heights and 
widths of the experimental data show good agreement 
with the limiting case theory when |l/fci?a| > 1. The 
width of the molecular peak is limited by the Fourier 
width of the Bragg pulse (~ 20 kHz). 

The integral of X{q,S) from equation ([T]) over 6 is 
proportional to the true static structure factor, S{q) = 
hN-^ J S{q,S)dS [H,!!^. Experimentally, if TBr is fixed 
and rJ^j. is known for a series of experiments, S{q) may 
be obtained to within a scaling factor by directly inte- 
grating a Bragg spectrum. 

In two component Fermi gases S{q) is determined by 
correlations between particles in all combinations of spin 
states. Labelling the two states with t and | arrows, 
S{q) = 5*11 (g) + S^iiq), wherein (g) = Sii{q) for 
N-^ = Ni and = [21j. Low energy exci- 

tations {q < kp) are suppressed by Pauli blocking and 
S-\i{q — * 0) — ^ 0. However, in our experiments q/kp ~ 5 
(> 1) and particles in the same state are essentially un- 
correlated so ^^^(g) 1. The total static structure fac- 
tor is then S{q) = 1 -|- S^i{q) which provides a direct link 
to correlations between spin-up/spin-down particles via 
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FIG. 3: Experimental static structure factor, Sexp(q = Skp) 
as a function of l/kpa and magnetic field (top). 
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the Fourier transform. For a uniform gas, 

(2) 

where n is the density and g^^ {r) represents the spin- 
up/spin-down two-body correlation function 21|. Quan- 
tum Monte Carlo simulations for a homogeneous gas pre- 
dict S{q > kp) to vary monotonically from 2 to 1 as 
1/kpa goes from being large and positive to large and 
negative In Fig [3] we show the experimentally deter- 
mined static structure factors, Sexp{q = 5kp), obtained 
by integrating the spectra in Fig. [^over 6. These are nor- 
malised so that Sexpiq = 5kp) at 750 G (l/kpa = 1.4) is 
two, corresponding to the bound molecule limit. S-\i{q) 
decays from the BEC to BCS sides of the Feshbach res- 
onance due to the decay of g^^'{r). 

For positive 1/kpa the measured Sexpiq = 5kp) ap- 
proaches the ideal molecular value of 2 closer to the Fes- 
hbach resonance than predicted for a uniform gas [2l| . 
Our experiments employ a nonuniform trapped gas and 
kp refers to the local Fermi wavevector at peak cloud 
density. Away from the trap centre kp{r) = [67r^n(r)]^/'^ 
decreases leading to a higher local interaction parameter 
l/kp{r)a. For a uniform gas, however, kp and 1/kpa 
are constant across the cloud. Applying equation [2] to a 
trapped gas using a local density approximation would 
involve averaging over regions with higher l/kp{r)a so a 
trapped gas would approach ideal molecular limit more 
rapidly than a uniform gas of equal 1/kpa. 

The existence of pairs at unitarity and above a Fesh- 
bach resonance is a quantum many body effect due to the 
strong interactions between the two species. As the den- 
sity is lowered, interactions decrease and pair correlations 
decay. To investigate this we apply a Bragg pulse reso- 
nant with pairs/molecules (5/27r = 145 kHz) to clouds 
which have been released from the trap and allowed to 
expand for a variable time. Imaging is always performed 
4 ms after the Bragg pulse and a reference image of an un- 
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FIG. 4: Response of trapped gases to a Bragg pulse resonant 
with pairs {5/211 = 145 kHz) as a function of cloud density. 
Magnetic field values and (l/kpa) are listed in the legend. 



scattered cloud is obtained for each pre-expansion time. 
The measured centre of mass displacements are plotted 
in Fig. m Peak densities uq were obtained from the ref- 
erence images and no x (67r^/fc|n)~^ can be larger than 
the ideal gas value of one because of interactions. 

On the BEC side of the Feshbach resonance the scat- 
tered pair signal starts high and decreases towards an 
asymptotic value as the density is lowered. During ex- 
pansion mean field energy is converted into kinetic energy 
and fewer molecules are resonant with the Bragg pulse 
[38| . A nonzero asymptotic value, when collisions are 
negligible, implies the scattered particles are real bound 
molecules. As the magnetic field is increased towards the 
Feshbach resonance, this asymptotic level decreases due 
to the stronger repulsive interactions. At unitarity, the 
scattered signal starts high and tapers off eventually dy- 
ing away to zero at densities below ~ 0.1(67r^/A;|,). This 
indicates that no bound molecules are present, however, 
in higher density clouds, pair correlations clearly exist. 
Similar behaviour is observed at 850 G and 890 G but the 
initial high density pair signal is lower. This technique 
allows us to distinguish between true bound molecules 
and paired atoms through the existence/nonexistence of 
a nonzero low density scattered signal. 

In summary, we have presented the first Bragg spectro- 
scopic study of a Fermi gas in the BEC-BCS crossover. 
Bragg spectroscopy allows a direct probe of two-body 
correlations and the transition from molecular to atomic 
behaviour has been observed. At unitarity and above the 
Feshbach resonance, pair scattering is observed. Pairing 
is seen to decay as the density is lowered, highlighting 
the many body nature of pair correlations. In future 
work, Bragg scattering could measure the temperature 
dependence of pairing and condensation 29]. The large 
momentum transfer used here prevents us from resolving 
the pairing gap. Howver, low (tunable) q Bragg scatter- 
ing could probe the superfluid excitation spectrum and 
the form of g^^ (r) throughout the crossover. 
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